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Abstract 

The formation of fibrillar aggregates seems to be a common characteristic of polypeptide chains, although the observation of 
these aggregates may depend on appropriate experimental conditions. Partially folded intermediates seem to have an important 
role in the generation of protein aggregates, and a mechanism for this fibril formation considers that these intermediates also 
correspond to metastable states with respect to the fibrillar ones. Here, using a coarse-grained (CG) off-lattice model, we carry 
out a comparative analysis of the thermodynamic aspects characterizing the folding transition with respect to the propensity 
for aggregation of four different systems: two isoforms of the amyloid /3-protein, the Src SIf3 domain, and the human prion 
proteins (hPrP). Microcanonical analysis of the data obtained from replica exchange method (REM) is conducted to evaluate 
the free-energy barrier and latent heat in these models. The simulations of the amyloid /3 isoforms and Src SH3 domain 
indicated that the folding process described by this CG model is related to a negative specific heat, a phenomenon that can 
only be verified in the microcanonical ensemble in first-order phase transitions. The CG simulation of the hPrP heteropolymer 
yielded a continuous folding transition. The absence of a free-energy barrier and latent heat favors the presence of partially 
unfolded conformations, and in this context, this thermodynamic aspect could explain the reason why the hPrP heteropolymer 
is more aggregation-prone than the other heteropolymers considered in this study. We introduced the hydrophobic radius of 
gyration as an order parameter and found that it can be used to obtain reliable information about the hydrophobic packing 
and the transition temperatures in the folding process. 

PACS numbers: 87.15. Aa,87.14.Ee,87.15.Cc,05. 70. Fh 
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I. INTRODUCTION 

Different diseases, including Alzheimer's [H [2] (AD), 
Parkinson's [3], and Huntington's [H [5] diseases, are 
known to be the result of neurodegenerative processes 
caused by formation of fibrillar aggregates of proteins, 
known as proteinopathies, which occurs primarily in 
the cellular structures [BUlOj. Because a number of 
proteins have been found to form fibrillar aggregates 
without sequence identity or structural homology under 
appropriate conditions in vitro, one can conclude that 
this is a common characteristic of polypeptide chains 
[3 Uni [H] , that wouU reflect a property of their main 
backbone. Therefore, the use of detailed interatomic 
potentials in computational experiments employed for 
the thermodynamic characterization of proteins with 
known aggregation propensity seems to be inessential and 
thus, some general features can be extracted with the 
aid of simpler models [12]. In particular, we followed ref. 
[13] and utilized a coarse-grained off-lattice model 14J in 
our simulations. In the simulations, the 20 naturally oc- 
curring amino acid residues are replaced with monomers 
whose hydrophobic-polar characters are denoted by A 
and B, respectively. Thus, our investigation of the ther- 
modynamic features of peptide chains with known amy- 



* Present address: Universidade Tecnologica Federal do Parana, 

Toledo, PR, Brazil; frigori@utfpr.edu. br^ 
t lerizzi@usp.br 
■t |alves@tt'clrp.usp.br| 



loidogenic propensity assumes that the nonbonded inter- 
actions are described by the hydrophobic-polar character 
of these monomers in all extension. The hydrophobicity 
of the side chains correlates with the aggregation rate 
[151 116] , so it is an important physicochemical ingredient 
promoting the nucleation of fibrillar aggregates. The mi- 
crocanonical analysis of the heteropolymer aggregation 
process has been performed [131 [T7] and supports the 
main idea that simple models encompassing hydrophobic 
interactions among monomers or segments present the 
overall aggregation phenomenon through a temperature- 
driven first-order transition. A more comprehensive 
study on the properties of the aggregation phenomenon 
shows that the change in the aggregation rate depends on 
diverse factors, grouped into extrinsic or intrinsic factors 
[T51 [T5] . Extrinsic factors include physicochemical prop- 
erties related to the polypeptide environment. The in- 
trinsic factors associated with amyloid formation refer to 
the characteristics of the polypeptide chains. Factors like 
hydrophobicity, propensity to form the a-helical struc- 
ture, propensity to form the /3-sheet structure, overall 
charge, and patterns of polar and nonpolar residues have 
been demonstrated to influence the aggregation rate. 
The propensity of globular proteins to adopt a different 
secondary structure content is another intrinsic factor 
related to the stability of the polypeptide chain. There 
is strong evidence that partially folded intermediates, 
such as for the A/3 peptide in the case of Alzheimer's, 
or with substantial unstructured conformations play an 
important role in the formation of amyloid fibril [191 - 
122]. Thus, conformational destabilization seems to be 



a natural requirement for polypeptide chains to achieve 
new arrangements even under physiological conditions, 
which culminates in their assembly of amyloid fibrils 
[53]. On the other hand, unstructured or destabilized 
conformations do not seem to be a necessary condition to 
promote aggregation [211 [23 , although results from mu- 
tations on the amyloid /3-protein (A/3) have shown that 
the less stable A/3 variants have the fastest nucleation 
kinetics in fibrillization [3H1 [57] . 

Here, thermodynamic features of the amyloid /3- 
peptides (^,940 and ^/342), the Src SH3 domain, and 
the human prion protein (hPrP) are compared in silico. 
Moreover, we try to associate those features with their 
aggregation propensity levels. The amyloid /3-protein is 
the principal component of the amyloid plaques found 
in AD patients [55]. It is a small peptide that is 
predominantly constituted by 40 or 42 residues. The 
Src SH3 domain in turn, is a globular domain with 56 
residues. Although this domain is non-pathogenic, it 
is expected to behave quite similarly to prion proteins 
during the aggregation transition. Actually, it has been 
experimentally shown that the SH3 domain aggregates 
to form amyloid fibrils [8] [29]. Molecular dynamics 
simulation of Src SH3 domains also led to aggregation 
conformations when the amino acid residue interactions 
were modeled by the Go potential |30j . The prion protein 
with PDB entry IHJM [31] is considered in this study. A 
large number of observations have led to the conclusion 
that partial unfolding of the monomeric cellular prion 
protein (PrP*^) seems to be necessary for the generation 
of the misfolded intermediates during fibrillization of the 
scrapie isoform PrP^'' [51 ] 155] , [55] . 

In addition, we emphasize that thermodynamic equi- 
librium properties are better analyzed in the microcanon- 
ical ensemble because it leads to the correct entropic 
characterization [34] and sets a simple way of determining 
free-energy barriers. In fact, many systems present 
equilibrium properties in the microcanonical ensemble 
but do not have their equivalents in the canonical en- 
semble [55H57] . This is a consequence of the noncon- 
cavity of the entropy function, which produces energy- 
dependent equilibrium states in the microcanonical en- 
semble that cannot be associated with any temperature- 
dependent equilibrium states in the canonical ensemble. 
These missed states are the nonequilibrium ones in the 
canonical ensemble, and they generally correspond to 
metastable or unstable states in this ensemble [38]. The 
most relevant feature that emerges from the nonconcavity 
of the entropy is the occurrence of a negative micro- 
canonical specific heat. Hence, it is preferable to perform 
a microcanonical analysis so that the metastable states 
that may arise as a consequence of the dynamic structural 
conversion between the native and unfolded states can be 
properly accounted for. In this work, we use the replica- 
exchange method (REM) [39^ and ST-WHAM-MUCA 
[40| analysis to obtain such microcanonical results and 
to base our discussions on. 



II. MODEL AND SIMULATION ALGORITHM 

To perform the simulations, the primary structure of a 
protein is mapped through the Roseman hydrophobicity 
scale [H] onto a string of elements A and B located at 
the Cq atoms. Table |l| includes the PDB codes used in 
this study and the corresponding AB sequences. The 
interaction for chains with N monomers is described by 
the following energy function [14] . 

N^2 N-2 N / , .\ 

k=l i=l j=i+2 \ ^3 ^3 I 

(1) 

Here, Ok is the angle between three successive monomers, 
and denotes the distance between the monomers i and 
j in the chain. The coupling constant C{ai,aj) in the 
Lennard- Jones type potential depends on the pairwise 
hydrophobic details, 

{1, fJi = aj ^ A 
0.5, a, = <Jj=B (2) 
-0.5, (7i ^ (Tj . 

The spherical-cap algorithm [42 was used in the simula- 
tions in order to update conformations. This algorithm 
utilizes spherical coordinates for the generation of new 
positions for each monomer in the chain. 

Conformations at different temperatures were obtained 
by REM. Figure [T] presents a log plot of the temperatures 
used in the replica exchange simulations for the het- 
eropolymers A/3, Src SH3, and hPrP. The temperature 
sets were determined by the following protocol: given 
an initial inverse temperature /3„, we determined the 
next inverse temperature Pn+i by simulating only two 
replicas of the system with the Metropolis algorithm and 
the AB force field. First of all, we performed REM using 
Ns sweeps and N swaps exchange conformation moves so 
as to equilibrate both replicas at the same reference 
temperature T„ = l//3„, with the Boltzmann constant 
ks — 1. Next, the inverse temperature of one replica 
was increased by a small variation S/3, and N^^aps x A^s 
updates were performed again trying to exchange replicas 
after Ng sweeps. If the fraction of accepted replica 
exchanges face was approximately equal to a probability 
Pacci this new inverse temperature Pn+i was accepted 
as reference; otherwise, the inverse temperature was 
increased by (5/3 once again. This procedure allowed 
for the recursive determination of all the temperatures. 
To compute facc^ ^swaps = 50 and Ng = 2000 were 
employed. The acceptance probability pacc was set to 
0.40 for the amyloid /3-peptides, and 0.30 for the Src SH3 
domain and the hPrP protein. It has been demonstrated 
that these Pacc values provide a convenient number of 
round trips between extremal temperatures [131 1^ . 
For the all systems, evaluation of the temperature sets 
started at the inverse temperature /3i — 0.5, with 
(5/3 — 0.01. A simple statistics with five independent 
simulations was used for face estimation, and thus for 
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establishment of the set of temperatures {Ti} to perform 
REM simulations for each heteropolymer. 

Final data production was obtained with 12 replicas 
for A/340 and A^42, 16 replicas for Src SH3, and 20 
replicas for hPrP. Simulations for data production were 
accomplished with Ng — 2000 sweeps and Ngwaps — 
10500 replica exchange moves, but for hPrP the number 
of exchange moves was doubled. Moreover, the above 
statistics was repeated five times, always starting from 
different initial conditions. Figure [2] shows conformations 
obtained from the AB model for the sequences in Table 
I, and sampled in the transition region. 

Data analysis was performed using the ST- WHAM 
procedure, which is an iteration-free approach that solves 
the WHAM equations in terms of intensive variables 
|45j . This procedure is a convenient way of obtaining 
the microcanonical inverse temperature estimates. On 
the basis of the assumption that the density of states can 
be obtained from 



H{E) 



(3) 



with H{E) being the energy histogram and W{E) the 
simulation weight, ST- WHAM states that the inverse 
temperature can be estimated from M independent sim- 
ulations {n — 1, • • • , M), 



dlnfliE) 
dE 



(4) 
(5) 



where /* = H^/Y^nHn, = dlnH^/dE, and /3f = 
—dlnWn/dE. The sum in n is over data produced 
by REM simulations, with the energy histograms 
and Wn = e~^/^" specified for each temperature r„. 
Thermodynamic quantities such as microcanonical en- 
tropy S{E) and specific heat Cy{E) are evaluated from 
a multicanonical entropy- like solution [ID], the so-called 
ST-WHAM-MUCA procedure. 



III. NUMERICAL SIMULATIONS AND RE- 
SULTS 

For comparative purposes, thermodynamic quantities 
such as the entropy, changes in the free energy A_F, 
and latent heat £ are given as a function of the specific 
energy e = E/N, where N stands for the total number 
of monomers in the system. These quantities are sum- 
marized in Table [III 

Because the force field treats the hydrophobic bonding 
interactions as the driving force in the folding process, it 
is important to analyze how the hydrophobic monomers 
behave as a function of the temperature. In fact, a 
recent study has shown that, besides hydrogen bonding, 
hydrophobicity is a fundamental interaction in the com- 
peting processes leading to folded or misfolded proteins 



[46) . As we will show, the folding behavior can also be 
investigated by taking into account the spatial distribu- 
tion of hydrophobic monomers by means of a radius of 
gyration restricted only to this type of monomer. It has 
been demonstrated that this hydrophobic radius rh is a 
suitable scoring function for discrimination between the 
native structure and other conformations [47[ . 



A. Microcanonical analysis of the A/3 models 

Figure |3]ja) exhibits the estimates of /3(e) for the 
A/340 heteropolymer. The microcanonical entropy S{e) 
is calculated using these ST- WHAM estimates of /3(e) 
according to the updating procedure ST-WHAM-MUCA. 
This entropy presents the so-called convex intruder (fig- 
ure not shown) and produces a negative specific heat in 
the energy range [—0.15,-0.07], which is a signature of 
a first-order phase transition. Figure depicts the 
behavior of Cv{e). A small positive peak in C,j(e) is 
observed at e = —0.579, presumably due to further com- 
paction of the heteropolymer chain. This latter transition 
is not related to a decrease in the microcanonical entropy 
with rising energy in this region. Therefore, /3(e) does not 
present a van der Waals-like curve in Fig. [sjjb), which 
results in a continuous phase transition at /3c = 2.695. 

To evaluate the free-energy profile as the heteropoly- 
mer assumes a folded conformation characterized by a 
stable phase within the metastable one, we considered 
the change of the microcanonical entropy between these 
phases in the vicinity of the transition temperature Tf — 
l/(3f. The energy region of interest is defined by the 
Maxwell construction in the caloric curve /3(e). This 
region is limited by the energies e^ and e?, {Sa < £b) where 
the horizontal straight line of the Maxwell construction 
intercepts the caloric curve /3(e) and identifies the inverse 
of the canonical first-order transition temperature /3y, 
which is related to the change in the entropy. 



1 S{eb) - S{ea) 

N Sb- £a 



(6) 



The energies Sa and £b define the discontinuity or the 
latent heat, observed in van der Waals loops. To calculate 
the change in the free energy AF(e) and /3/ from S'(e), 
we defined a shifted entropy between £a and e^, AS'(e) — 
N{A -f /3/e) — '5'(e), where A and /3/ are such that a 
canonical entropy, Scani^) — N{A + /3/e), is defined in 
that energy range. The constants A and Pf are such 
that Scan{£a) = S{ea) and Scan{£b) = S{£b)- These 
conditions yield £a = —0.191 and e^ — —0.026, at 
Tf = 0.692(1) for A/?40. 

The corresponding change in the free energy can be 
estimated as 



f3f[F{£)-Fa]^f3fN£-S{£), 



(7) 



where /3/ = /?/ N£a — S{£a) is the reference free energy. 
Our results for AF(e) are shown in the inset of Fig. 
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[sjja). This analysis yields the free-energy barrier AF = 
0.038(2) at the folding temperature Tf. Latent heat is 
a consequence of the free-energy barrier that prevents 
the system from moving to a stable conformation in the 
new phase. Therefore, the smaller the latent heat, the 
higher the probability that a spontaneous fluctuation 
will give rise to this new phase. The estimate for the 
latent heat per monomer associated with this transition 
is ^ = 0.165(3). 

A similar analysis follows for the heteropolymer ob- 
tained with the PDB entry IZOQ for the Af342 peptide. 
Results are presented in Fig. [3] However, for this het- 
eropolymer, the caloric curve displays a less pronounced 
van der Waals-like behavior around e = —0.2 as a conse- 
quence of a minor nonconcavity of the microcanonical 
entropy. This model yields AF = 0.014(1) (inset of 
Fig. Qa)) and latent heat £ = 0.125(2). Therefore, the 
free-energy barrier separating the folded and unfolded 
states of the heteropolymer describing A/342 is smaller 
compared with the A/340 heteropolymer. This corre- 
sponds to a less severe restriction to possible movements 
returning the heteropolymer to intermediate conforma- 
tions, depending on how stable the folded conformations 
are |3S]. A small positive peak in the specific heat at 
e = —0.686 signals a continuous transition (/3c = 2.632) 
similar to the one observed for the A/340 heteropolymer. 

The results of the hydrophobic gyration radius for 
the A/340 and A/342 heteropolymers as a function of 
temperature are given in Fig. [5] The folding transition is 
expected to be accomplished by fast change in the spatial 
distribution of the monomers. In particular, this change 
is expected to be highly sensitive to the hydrophobic 
monomers, as illustrated in Fig. [sja). Thus, we 
hypothesized that the temperature where the derivative 
< drh/dT > hits its maximum represents the occurrence 
of a thermodynamic transition. These temperatures, 
denoted by Tr, are easily identified in Fig. [sjjb): Tr = 
0.68(1) and = 0.76(1) for A/340 and A;942, respec- 
tively. These figures also identify the latter compaction 
transitions at temperatures 0.37(1) and 0.38(1) for the 
A/340 and A/342, respectively. Estimates of Tr are in very 
good agreement with the values Tf obtained via Maxwell 
construction for both heteropolymers. This indicates 
that rh can be considered a reliable order parameter 
for depiction of a hydrophobic profile as a function of 
temperature. This conclusion is also supported by results 
attained for the Src SH3 and hPrP heteropolymers as 
shown in the following subsections. 



B. Microcanonical analysis of the Src SH3 model 

The SH3 domains have attracted much interest be- 
cause they represent typical examples of proteins that 
fold via a two-state mechanism. It is largely accept that 
the physical process underlying protein folding in these 
cases is based on the nucleation-condensation scenario 
[3nHS2]. In fact, a sharp transition at the folding 



temperature Tf between the unfolded and folded states 
has been observed [49l [50] . 

Results from simulations for the Src SH3 heteropoly- 
mer obtained with the AB model are shown in Fig. 
|6] As expected for such peptide with a clear two-state 
transition mechanism, a van der Waals-like loop can be 
detected in the caloric curve (Fig. [6]ja)). Therefore, the 
microcanonical specific heat presents the typical behavior 
observed in a first-order phase transition (Fig. [6]^b)). 
The Maxwell construction leads to Sa — —0.182 and 
Eb = 0.0079, with Tf = 0.658(1). For this model, we 
obtained a free-energy barrier AF = 0.068(2) (inset in 
Fig. |6][a)) and a latent heat per monomer i = 0.190(8). 
Two small positive peaks are seen in C„(e) at e = —0.21 
and £ = —0.70. Both transitions are likely related to 
compactions. The behaviour of and its derivative 
are included in Fig. ^ The maximum of < drh/dT > 
occurs at Tr = 0.66(1). This feature also signals the 
existence of a second (minor) peak in the derivative of 
at Tr = 0.47(3), which is related to the transition 
observed in Fig. |6ja) at e = —0.70. 



C. Microcanonical analysis of the hPrP model 

The next heteropolymer corresponds to the extracellu- 
lar globular domain hPrP with 104 residues. Compared 
to the previous models, the numerical results do not 
furnish any convex intruder in S'(e). The behavior of 
(3{e), displayed in Fig. [Tj^a), does not produce a van 
der Waals-like curve. The microcanonical analysis yields 
positive specific heats associated with continuous phase 
transitions (Fig. [7][b)). The maxima of Cy{e) occur 
at El = —0.36 and £2 = —0.064, which correspond to 
the phase transition temperatures Ti = 0.596(2) and 
T2 = 0.697(1), respectively. The behaviour of and 
its derivative are also included in Fig. [5j Here, there is a 
clear transition at Tr — 0.60(1) and just a tiny change in 
the derivative of the hydrophobic radius at Tr = 0.32(1). 

Since our findings for this heteropolymer did not 
evidence a reduction in the microcanonical entropy as 
we moved toward the unfolded states, we may argue 
the following. Either the free-energy barrier separating 
the folded and unfolded states for this heteropolymer is 
not large enough to be revealed by our potential energy 
function because it does not contain interactions that 
reproduce the native contacts, or in fact the transition 
does not occur via a two-state mechanism. Interest- 
ingly, contradictory observations in what concerns the 
prion protein folding have been reported from diverse 
experiments. Experimental observations indicate that 
the native folding pathway involves the two-state process 
despite other evidences to the contrary (see, for example 
(33j and references therein). 
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IV. DISCUSSION AND CONCLUSIONS 

The ability of peptide chains to form fibriUar ag- 
gregates seems to foUow from their propensity to ag- 
gregate under conditions that permit partial unfolding. 
These fibrillar structures originate from misfolded protein 
chains via a complex process. A mechanism for the 
formation of the fibrils considers that the intermediates 
correspond to metastable states with respect to the 
fibrillar states. In this sense, the existence of these 
aggregation-prone states can be a consequence of the 
low degree of stability of the native states. It has 
been demonstrated with a simple lattice model [55] 
that low-energy conformations populate the aggregation- 
prone states. Such conformations are identified in the 
variety of lowest energy oligomers and protofilaments 
obtained when multiple chains are present in the three- 
dimensional lattice. Here, we performed REM simu- 
lations of four biologically inspired heteropolymers to 
relate the thermodynamic properties with their aggrega- 
tion propensities. Considering the amyloid /3-peptides, 
there is experimental evidence that the A/342 is more 
aggregation-prone than A/340 |53j . Our findings for 
the free-energy barriers and latent heats of these het- 
eropolymer systems, as listed in Table [llj indicate a two- 
state folding process. Smaller AF and £ for the A/342 
heteropolymer, compared with the respective values for 
A/340, indicate that the formation of native states are 
facilitated for A/342. On the other hand, these numerical 
values signal a weaker first-order transition for A/342. 
This may give higher chances for the A/342 native-like 
conformations to adopt partially folded intermediates 
under similar stability conditions. As a matter of fact, 
monomeric A^ peptides are considered to be intrinsically 
disordered and therefore, A/342 can cross the free-energy 
barrier to misfold conformations more easily |26|. 

For the Src SH3 heteropolymer, our results show 
a stronger first-order phase transition, confirming the 
two-state character expected for this system. Thus, 
on comparative grounds, it is reasonable to understand 
the experimental requirements to destabilise the native 
conformation [54] to produce some partially unfolded 
conformations as a prerequisite for its self-assembly. 

The confiicting evidence about the mechanism of the 
folding process involving PrP shows how difficult the 
experimental measurements are [33J. Interestingly, our 
in silico experiment suggests that if there is a free- 



energy barrier separating the folded and unfolded states 
of the hPrP heteropolymer, it is not large enough to be 
revealed by our simple potential energy function. We 
argue that the absence of a free-energy barrier favors 
the presence of partially unfolded conformations in the 
hPrP heteropolymer, which could explain why it is more 
aggregation-prone than the other heteropolymers. Of 
course, our conclusions are based on a simple force field 
which was not designed, for example, to reproduce the 
native state. However, we expect that the hydrophobic 
force incorporates the main aspects that produce the 
transition-state configurations. 

Because the propensity of peptides and proteins to 
form aggregates depends greatly on the sequence, Pawar 
et al. |18j defined a phenomenological equation to express 
this feature. We calculated the intrinsic score Zagg 
of aggregation for our polypeptide chains at pH 7 by 
means of the Zyggregator algorithm [56' . This algorithm 
yielded the values presented in Table [TTl Observation of 
the Z-scores reveals that the hPrP polypeptide is more 
aggregation-prone than the other polypeptide chains. 
This conclusion agrees with the thermodynamic phase 
transition results listed in Table Hi] for these different 
heteropolymer systems, if we consider that weaker phase 
transitions facilitate the coexistence of mixed conforma- 
tions. 

We also computed the hydrophobic radius of gyration 
and as can be seen from Fig. [5] it is a convenient 
order parameter for identification of the folding tempera- 
tures. The agreement with the microcanonical estimates 
demonstrates that the spatial packing of hydrophobic 
monomers furnishes reliable information about the fold- 
ing process. More importantly, this quantity can easily 
be analyzed in any protein folding study and does not 
depend on any other information about the protein 
chain like the usual reaction coordinates needing some 
information about native contacts. 
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FIG. 1. Temperatures used to perform REM simulations. 




FIG. 2. Illustrative sampling of conformations for the AB sequences representing (a) A/?40, (b) A/342, (c) Src SH3, and (d) 
hPrP proteins in the transition region. Dark (in color: red) color indicates hydrophobic monomers. (These images were made 
with VMD software support). 
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FIG. 3. (a) Microcanonical inverse temperature estimates /?(e) = 1/T{e), and (b) microcanonical specific heat C„(e) = 
—0^/{dp/de) for the A/340 heteropolymer. The inset of figure (a) shows the free-energy changes constructed from S{e) at j3f. 




1.5 h 
1 

24 
12 


-12 
-24 







-(b) 

. I.I.I. 


1 . 1 



-1 -0.8 -0.6 -0.4 -0.2 

e 



FIG. 4. (a) Microcanonical inverse temperature estimates /3(e), and (b) microcanonical specific heat C„(e) for the A/342 
heteropolymer. The inset of figure (a) shows the free-energy changes constructed from S{s) at /3/. 
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FIG. 6. (a) Microcanonical inverse temperature estimates /5(e), and (b) microcanonical specific heat C„(£) for the Src SH3 
heteropolymer. The inset of figure (a) shows the free-energy changes constructed from S{e) at /3/. 



9 




-1 -0.8 -0.6 -0.4 -0.2 
8 



FIG. 7. (a) Microcanonical inverse temperature estimates /3(e), and (b) microcanonical specific iieat Cv(e) for tfie fiPrP 
lieteropolymer. 



TABLE 1. PDB codes and Roseman mapping. 



Protein PDB code Residues Sequence 



amyloid f5 


2LFM 


40 


BABAB 


BBBBB 


BABBB 


BAAAA 


ABBAB 


BBBBA 


AABAA 


ABBAA 






amyloid /3 


IZOQ 


42 


BABAB 


BBBBB 


BABBB 


BAAAA 


ABBAB 


BBBBA 


AABAA 


ABBAA 


AA 




Src SH3 


INLO 


56 


AAAAA 
BBAAA 


BBBBB 
B 


BABAB 


ABABB 


BBBAB 


AABBA 


BBBAA 


AABBA 


AABBA 


BBAAB 


prion protein 


IHJM 


104 


ABBBA 


ABBAA 


BBAAA 


BABBB 


BBBBB 


BBBBA 


BBBAB 


BABBB 


AABBB 


BBBBB 








AABBA 


ABAAA 


BBBAA 


AAAAB 


BBBAA 


BABAB 


AABBA 


ABBAA 


AABBB 


BBBBA 








BBBB 





















TABLE II. Comparative results for the heteropolymer models. 



Model 


Tf 


AF 




T, 


Zagg score [56] 


A/?40 


0.692(1) 


0.038(2) 


0.165(3) 


0.68(1) 


0.90 


A/?42 


0.769(1) 


0.014(1) 


0.125(2) 


0.76(1) 


0.94 


Src SH3 


0.658(1) 


0.068(2) 


0.190(8) 


0.66(1) 


0.96 


hPrP 


0.596(2) 






0.60(1) 


1.14 
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